Massive graviton as a testable cold dark matter candidate 
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We construct a consistent model of gravity where the tensor graviton mode is massive, while 
linearized equations for scalar and vector metric perturbations are not modified. The Friedmann 
equation acquires an extra dark-energy component leading to accelerated expansion. The mass of 
the graviton can be as large as ~ (10 15 cm) -1 , being constrained by the pulsar timing measurements. 
We argue that non-relativistic gravitational waves can comprise the cold dark matter and may be 
detected by the future gravitational wave searches. 
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1. Introduction. The current cosmological model is in 
a beautiful agreement with the data [1]. However, it re- 
quires introduction of exotic density components (dark 
matter, dark energy) with abundances highly tuned to 
the baryonic matter. This motivates interest in modified 
theories of gravity deviating from the Einstein theory at 
large distance scales. Generically, in such theories the 
graviton has a non-zero mass. The common lore is that 
the inverse graviton masses significantly smaller than the 
current Hubble scale are not phenomenologically allowed. 
In this paper we demonstrate that the inverse graviton 
mass can be not only significantly smaller than the cur- 
rent size of the Universe, but even many orders of mag- 
nitude smaller than the galactic scales. We argue that 
massive graviton provides specific signatures for gravita- 
tional wave experiments and may even account for the 
cold dark matter (CDM) in the Universe. 

Recent studies of the Fierz-Pauli theory of massive 
gravity [2] and brane world scenarios where the four- 
dimensional graviton has a non-zero mass [3, 4] strongly 
suggest [5-11] that Lorentz-invariant models of massive 
gravity suffer either from the presence of ghosts (fields 
with a wrong sign of the kinetic term) , or from the vDVZ 
discontinuity due to extra graviton polarizations [12, 13] 
and strong coupling at the low energy scale. It is possible 
that the account for the effects of local curvature may 
solve these problems in some models [14-17]. Another 
possibility which attracted attention very recently [18- 
22] is to allow for a violation of Lorentz invariance. In 
particular, a class of models was found [22] where tensor 
graviton mode is massive, vDVZ discontinuity and strong 
coupling problems are absent, while the absence of ghosts 
and rapid classical instabilities is ensured by the residual 
reparametrization symmetry 
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x l being the spatial coordinates. These models are the 
focus of the current paper. 

2. The model. In the covariant formalism of Ref. [22] 
(see also Refs. [5, 19]), the action for the theory of massive 
gravity contains the metric g^ v and four scalar Goldstone 



fields <fP , 4> l (i = 1, . . . , 3). In the presence of the residual 
symmetry (1) it reads 

S = J d 4 x^ [~M 2 Pl R + k 4 F(X, W ij , ...)], (2) 

where X and W % i are the scalar quantities constructed 
from the Goldstone fields and the metric, 
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and F is a function to be constrained later. We assume 
that the Goldstone sector is characterized by a single en- 
ergy scale A. Dots in Eq. (2) stand for higher-derivative 
terms. Latin indices i,j are contracted using 5y. 

We require that the model admits a background solu- 
tion with the metric g^ equal to the Minkowski metric 
r\y, v and the scalar fields taking the form 
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for some constants a and b. For a generic function F 
such a solution always exists. In the "unitary gauge" 
where the Goldstone fields are fixed to their vacuum val- 
ues (4), the second term in the action (2) gives rise to 
the following mass term for the metric perturbation h^ Vl 
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where the values of the mass parameters m a are deter- 
mined by the first and the second derivatives of the func- 
tion F(X, W lj ) at the vacuum values of its arguments as 
defined by eqs. (3) and (4). The overall scale m of the 
graviton masses is related to A as m ~ k 2 /Mpi. The 
analysis of Ref. [22] implies that A plays the role of the 
cutoff scale of the theory with the action (2). 

The residual reparametrization symmetry (1) arises in 
the unitary gauge as a consequence of the global sym- 
metry (j) 1 — > (jf + £, l (4>°) of the covariant action (2). This 
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symmetry implies, in particular, that there is no graviton 
mass term proportional to h^. 

As is usual in the linearized theory, it is convenient 
to consider separately tensor, vector and scalar metric 
perturbations (cf. Refs. [20, 22]). The tensor modes - 
transverse traceless gravitational waves hj? — have non- 
zero mass equal to to 2 [20] . There arc no propagating de- 
grees of freedom in the vector sector [22]. Moreover, the 
contribution of the mass term (5) in the vector sector has 
the form of a gauge fixing. Consequently, no modifica- 
tion of gravity arises in the vector sector at the order we 
are working. Finally, the energy-momentum tensor 8T^ V 
induces the following perturbations in the scalar sector, 

* = *b, (6) 
m 2 [3mj-m 2 (3m 2 -m 2 2 )} 1 <5r 00 
9E+ mi-ml{ml-ml) df Mf, t ' ^ 

where "J, $ are the gauge- invariant scalar potentials de- 
fined in a standard way [23], and ^e, &e are their values 
in the Einstein theory. The modification of gravity man- 
ifests itself in the last term in Eq. (7) . There is no vDVZ 
discontinuity as this term vanishes in the limit when all 
graviton masses uniformly go to zero. 

The extra term in Eq. (7) grows linearly with the dis- 
tance from the source, indicating the breakdown of the 
linearized theory. This growth cannot be eliminated by 
a proper choice of the gauge as $ is the gauge-invariant 
quantity. However, the Riemann curvature associated 
with the extra term goes to zero as 1/r at large r, so the 
space-time becomes flat far from the source. (This break- 
down of perturbation theory is very different in nature 
from the seemingly similar problem in the Fierz-Pauli 
theory [14], where it happens in the vicinity of the source. 
The close analogue of the phenomenon discussed here is 
the breakdown of perturbation theory far from the source 
in the three-dimensional classical Yang-Mills theory.) In 
the region where the non-standard term in Eq. (7) is still 
small it produces the r-independent force, imitating the 
effect of a halo with the density profile oc r _1 . 

The analysis of Eq. (7) in the region where it enters the 
non-linear regime goes beyond the scope of this paper. 
Instead, we chose the masses m a in such a way that the 
second term in Eq. (7) vanishes. It is important that 
this can be achieved by imposing, in addition to (1), the 
following dilatation symmetry, 

t -» At , x l -> \^x\ (8) 

where 7 is a real constant. At the linearized level this 
symmetry implies the following relations among masses, 

TOq = — 37TO4 , m\ — 3TO3 = 7~ 1 77i4, (9) 

which lead to the cancellation of the second term in 
Eq. (7) for any 7. Thus, when the symmetry (8) is im- 
posed, the only modification of gravity at the linearized 
level is the non-zero mass of the graviton. 

The inclusion of higher-derivative terms in the action 
(2) leads in general to the appearance of the dynamical 



degree of freedom in the scalar sector [22] . This degree of 
freedom is similar to that present in the ghost condensate 
model [19]. It has a healthy kinetic term provided the 
following inequality holds [22] , 
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The latter condition is compatible with Eqs. (9) and 
the requirement that the graviton mass is not tachy- 
onic, m\ > 0. The effects related to this degree of 
freedom are characterized by the huge retardation time 
~ m~ 1 (Mpi/ 'A) [19, 24, 25]. This time is larger than the 
current age of the Universe for the values of the graviton 
mass m specified below, so we can consistently neglect 
these effects. 

In the covariant formalism the residual symmetry (8) 
translates into the following global symmetry of the Gold- 
stone sector, <p° — > \<fi° , <j) % — > A _7 </>\ The action in- 
variant under the symmetries (1), (8) has the form (2) 
with the function F depending on the single combination 
X n W % K The case of the ghost condensate [19] emerges 
in the limit 7 — > 00 and requires a fine-tuning of F to 
obtain the Minkowski vacuum. The Minkowski vacuum 
with the scalar vev's of the form (4) exists for a general 
function F if 7 = 1 /d, where d = 3 is a number of spatial 
dimensions. For definitcness, in what follows we consider 
the case F = F(X 1 /' i W io ). 

3. Cosmological solutions. The spatially flat homoge- 
neous cosmological ansatz is 

ds 2 = a 2 (r])(dn 2 - dxfj , (11) 
0° = 0(7?), f =AV. (12) 

In what follows we assume that the rate of the expan- 
sion is much smaller than the energy scale A, so one can 
neglect higher derivative terms in the action (2). For sim- 
plicity, let us also assume that the function F depends 
only on the combination Z = X 1 ^ d 'W l Hij. The Einstein 
equations are reduced to the Friedmann equation 

{i) 2 -jkM'- + ^ nz)z - k ' nz) )- <13) 

where p m is the energy density of matter, and the field 
equation for <f)°, 

d n (a 3 F , (Z)WX~ 1 ^ = 0. (14) 

Eq. (14) implies Z = const or, equivalently, cx 
J di]a 4 (n). Then Eq. (13) takes the form of the standard 
Friedmann equation with the value of the cosmological 
constant determined by the value of Z, i.e. by the initial 
conditions in the Goldstone sector. Note that these initial 
conditions may be different in different regions of space. 
Therefore, this model is an example of the setup where 
de Sitter solutions with different expansion rates exist for 
any value of the vacuum energy. This property is a wel- 
come feature for the application of the weak anthropic 
principle [26] to the cosmological constant problem. 
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To summarize, we have constructed a consistent model 
where gravitational waves are massive, while linearized 
equations for the metric perturbations in the scalar and 
vector sectors, as well as spatially flat cosmological so- 
lutions, are the same as in the Einstein theory. In this 
model, the tests of (linear) gravity based on the solar 
system and Cavendish-type experiments [27] are auto- 
matically satisfied, while the main constraints are coming 
from emission and/or propagation of gravitational waves. 

4- Relic gravitational waves. Observations of the slow 
down of the orbital motion in binary pulsar systems [28] 
imply that the mass of the gravitational waves cannot be 
larger than the frequency of the waves emitted by these 
systems. The latter is determined by the period of the 
orbital motion which is of order 10 hours, implying the 
following limit on the graviton mass, 
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Let us estimate the cosmological abundance of relic 
gravitons. For this purpose we consider the transverse 
traceless perturbation of the metric h^. The quadratic 
action for hij in the expanding Universe takes the follow- 
ing form, 

M 2 Pl J d'kdrja 2 ^) (hi - (dkhijf - m 2 2 a 2 (ri)h 2 j) . 

(16) 

This has a form of the action for a minimally coupled 
massive scalar field. Therefore, gravitons in our model 
are produced efficiently during inflation (cf. Ref. [29]). 

To be concrete, consider a scenario where the Hubble 
parameter Hi is constant during inflation. This scenario 
may be realized, for instance, in hybrid models of infla- 
tion [30]. First, we need to check that the phenomeno- 
logically relevant values of parameters correspond to the 
regime below the cutoff scale of the effective theory, i.e. 
Hi < A. For the energy scale of inflation Ei ~ y/HiMpi 
this implies 



Ei < m^M%* 
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This value is high enough to allow for a successful baryo- 
genesis even for graviton masses of the order of the cur- 
rent Hubble scale. 

Consider now the production of massive gravitons. As- 
suming the above scenario of inflation, the perturba- 
tion spectrum for the massive gravitons is that for the 
minimally coupled massive scalar field in the de Sitter 
space [31], 
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Supcrhorizon metric fluctuations remain frozen until the 
Hubble factor becomes smaller than the graviton mass, 
when they start to oscillate with the amplitude decreas- 
ing as a~ 3 / 2 . The energy density in massive gravitons at 
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where we integrated in Eq. (18) over the modes longer 
than the horizon. Today the fraction of the energy den- 
sity in the massive gravitational waves is 
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where z a is the redshift at the start of oscillations, 
H <~ m 2 is the Hubble parameter at that time, H e w 
0.4 • 10~ 12 s _1 is the Hubble parameter at the mat- 
ter/radiation equality, and z e w 3200 is the correspond- 
ing redshift. Combining all the factors together one gets 



O s ~ 3 • 10 3 (m 2 • 10 15 cm) 1/2 I — 



(21) 



This estimate assumes that the number of c-foldings dur- 
ing inflation is large, In N e > H 2 /m 2 , which is quite nat- 
ural in the model of inflation considered here. 

According to Eq. (21), the massive gravitons are pro- 
duced efficiently enough to comprise all of the cold dark 
matter, provided the value of the Hubble parameter dur- 
ing inflation is about one order of magnitude below the 
scale A. We find it encouraging that one obtains fl g <~ 1 
when the initial energy density in the metric perturba- 
tions is close to the cutoff scale, pV^ ~ A. This suggests 
that other mechanisms of production unrelated to infla- 
tion (e.g., similar to those invoked for the axion or Polony 
fields) may naturally lead to the same result, Cl g <~ 1. 

The produced gravitons may cluster in galaxies. To 
account for the dark matter in galactic halos the graviton 
mass should satisfy (m») _1 < 1 kpc <~ 3 • 10 21 cm, where 
v ~ 10~ 3 is a typical velocity in the halo. 

5. Detection. Let us now briefly describe potential ob- 
servational signatures of the above scenario. Note first 
that at distances shorter than the wavelength, the effect 
of a transverse traceless gravitational wave on test mas- 
sive particles in Newtonian approximation is described 
by the acceleration hijX^ /2 (see, e.g., Ref. [32] for a re- 
view). The same is true for massive gravitational waves, 
the only difference being that the wavelengths are longer 
in the non-relativistic case, so the Newtonian description 
works for the larger range of distances. Thus, the non- 
relativistic waves act on the detector in the same way as 
massless waves of the same frequency. 

Let us estimate the amplitude of the gravitational 
waves assuming that they comprise all of the dark mat- 
ter in the halo of our Galaxy. The energy density in 
non-relativistic gravitational waves is of order M Pl rn\h 2 ^ . 
Equating this to the local halo density one gets 
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At the frequencies 10~ 6 H-10~ 5 Hz this value is well above 
the expected sensitivity of the LISA detector [33]. Note 
that in the close frequency range 1CP 9 -j- 10~ 7 Hz there 
is a restrictive bound [34] at the level il g < 1CT 9 on the 
stochastic background of the gravitational waves coming 
from the timing of the millisecond pulsars [35]. So, it is 
possible that our scenario can be tested by the reanalysis 
of the already existing data on the pulsar timing. 

The relic abundance of gravitons may depend on both 
the specific inflationary model and the details of the (un- 
known) UV completion of massive gravity. In general, 
massive gravitons may not comprise the whole of the 
CDM in the galaxy halos. It is important that the ex- 
pected LISA sensitivity allows to detect the presence of 
massive gravitons at the significantly lower level than in 
Eq. (22). 

6. Concluding remarks. In this paper we limited our- 
selves to a specific choice of the parameters (graviton 
masses and the constant 7 entering Eq. (8)) such that 
there is no modification of the Newton potential at the 
linear level, and the cosmological evolution remains stan- 
dard. We also did not consider possible non-linear effects, 
which may become a necessity at different choice of the 



parameters. A number of interesting questions is related 
to these effects, including the limits on graviton masses, 
clustering of massive gravitons in haloes and proper mod- 
ifications of Eqs. (13) and (14) to account for the direct 
coupling between Goldstone fields and gravitons. We ex- 
pect, however, that our main conclusions — that gravi- 
tons may have large masses and may be produced with 
cosmologically significant abundance — are generic in 
this class of models. In the relevant range of parame- 
ters, a specific signature of the gravitons with non-zero 
mass is a strong monochromatic signal in the detectors 
of gravitational waves. An independent measurement of 
the graviton mass may be performed at future gravita- 
tional wave detectors (for a review, see, e.g. [36]) operat- 
ing at higher frequencies by testing the delay between the 
electromagnetic and gravitational signals from a distant 
supernova explosion. 
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